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Abstract 

We consider the first LHC data for pp collisions in the framework of Regge phe- 
nomenology, and the Quark-Gluon String Model, and we present the corresponding 
predictions for both the integral cross sections and the inclusive densities of secon- 
daries, that are determined by Pomeron exchange. All parameters were fixed long time 
ago for the description of the data at fixed target energies. The first measurements of 
the total inelastic cross section by ATLAS and CMS Collaborations are in agreement 
with our calculations. Also the inclusive densities measured in the central region are 
in agreement with our theoretical predictions on the accuracy level of about 10—15%. 
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1 Introduction 



In Regge theory the Pomeron exchange dominates the high energy soft hadron inter- 
action, the contributions of all other exchanges to the total or inelastic cross sections 
becoming negligibly small at LHC energies. 

The Quark-Gluon String Model (QGSM) [Ij is based on Dual Topological Uni- 
tarization (DTU), Regge phenomenology, and nonperturbative notions of QCD. This 
model is successfully used for the description of multiparticle production processes in 
hadron-hadron [21 El IH El El 13 E], hadron-nucleus [H HU], and nucleus- nucleus [TT] 
collisions. In particular, the inclusive densities of different secondaries produced in pp 
collisions at ^/s = 200 GeV in midrapidity region were reasonably described in ref. [6] 
(see also refs. [Zl E])- 

In the QGSM high energy interactions are considered as proceeding via the exchange 
of one or several Pomerons, and all elastic and inelastic processes result from cutting 
through or between Pomerons [12]. Inclusive spectra of hadrons are related to the 
corresponding fragmentation functions of quarks and diquarks, which are constructed 
using the Reggeon counting rules [13]. The quantitative predictions of the QGSM 
depend on the values of several parameters that were fixed at once by the comparison 
of the corresponding calculations to the experimental data obtained at fixed target 
energies. 

The experimental data obtained at LHC allow us to test the stability of the QGSM 
predictions and of the values of the model parameters. Fortunately, one can see that 
the model predictions are in reasonable agreement with the first LHC data. 

In this paper we discuss the cross sections, inclusive densities of secondaries, and 
the antibaryon/baryon ratios at LHC energies in the framework of QGSM. 

2 Cross sections 

For the Pomeron trajectory 



where 7 = gi{0) ■ 02(0) is the Pomeron coupling, sq — 1 GeV^, and ap rises with energy 
as s^. The correct asymptotic behavior aj^^ ~ In^ s, compatible with the Froissart 



ap{t) = 1 + A + a'p-t , A>0, 
the one-Pomeron contribution to aj^^ equals 

(Tp = Svr ■ 7 ■ e^'^, with ^ = In s/sq , 



(1) 



(2) 
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bound [H], can be obtained by taking into account the multipomeron exchanges in the 
diagrams of Fig. 1. 




Figure 1: Regge-pole theory diagrams: (a) single, (b) double, (c) and triple Pomeron exchange in 
elastic hN scattering. The hadron couplings to two and three Pomerons are shown by dashed boxes. 

A simple quasi-eikonal treatment [15] yields to 

Of 
C 
1 

k ■ k\ ^ ' z Jo X 



<N = ■ f{z/2) , at^ = ^ • [f{z/2) - f{z)] , (3) 

/(^) = £ At ■ = 7 r ? ■ (1 - ^"^) ' (4) 



k=l 

2C-7 



z 



■e^^ , \ = R^ + a'p-i. (5) 

A 

Here, is the radius of the Pomeron and C is the quasi-eikonal enhancement coeffi- 
cient. 

At asymptotically high energies (2; ^ 1) we obtain 

tot _ 87i-a'p-A 2 .1 _ 47r ■ Qp ■ A . ff^^ 

'-'hN — Q ■ S ; '^hN ~ (J2 ' ^ ^ ^ ' 

according to the Froissart limit |14] . 

The numerical values of the Pomeron parameters were taken [3J to be : 
A = 0.139, a^ = 0.21GeV-^ 7 = 1.77GeV-^ i?^ _ 3.18 GeV-^ C = 1.5. (7) 

The predictions of Regge theory obtained wity these values of the parameters (and 
accounting for a small contribution from non-Pomeron exchange) are presented in Ta- 
ble. 1. 
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a*"^' (ATLAS [16j) 


(T*"^'(CMS [17J) 


900 GeV 


67.4 


13.2 


54.2 






7 TeV 


94.5 


21.1 


73.4 


69.4 ±2.4 


66.8 - 74.8 


14 TeV 


105.7 


24.2 


81.5 







Table 1. The Regge theory predictions for total, total elastic, and total inelastic cross sections (in 
mb) in pp collisions at LHC energies. 

These results are in agreement with those in ref. [TSj . 

The experimental points for a^^^ measured by ATLAS and CMS Collaboration [Ml 
[T7] presented in Table 1 (we omit the error bar coming from extrapolation) are in a 
reasonable agreement with our calculations. 

However, in the complete Reggeon diagram technique [19] not only Regge-poles 
and cuts but also more complicated diagrams (e.g. the so-called enhanced diagrams) 
should be taken into account. In the numerical calculations of such diagrams some 
new uncertainties appear due to the fact that the vertices of the coupling of n and m 
Reggeons are unknown. 

The numerical calculations which account for enhanced diagrams [201 1211 122] lead 
the values of a^^^ of the same order (± 10% at ^/s = 14 TeV) as those presented in 
Table 1. The values of a^^*^^ calculated in refs. [201 EI] are slightly smaller than the 
experimental values [T6l IT7] , and while two sets of calculations in ref. [22] lead to too 
large cross section, the third one [22] practically coincides with our results. 

3 Inclusive densities 

The Quark-Gluon String Model (QGSM) [H 121 [3] allows us to make quantitative pre- 
dictions of different features of multiparticle production, in particular, the inclusive 
densities of different secondaries both in the central and in the beam fragmentation 
regions. In QGSM high energy hadron-nucleon collisions are considered as taking place 
via the exchange of one or several Pomerons, all elastic and inelastic processes resulting 
from cutting through or between Pomerons |12j . 

Each Pomeron corresponds to a cylindrical diagram (see Fig. 2a), and thus, when 
cutting one Pomeron, two showers of secondaries are produced as it is shown in Fig. 2b. 
The inclusive spectrum of a secondary hadron h is then determined by the convolution 
of the diquark, valence quark, and sea quark distributions u{x,n) in the incident parti- 
cles, with the fragmentation functions G^{z) of quarks and diquarks into the secondary 
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hadron h. These distributions, as well as the fragmentation functions, are constructed 
by using the Reggeon counting rules [13]. Both the diquark and the quark distribution 
functions depend on the number n of cut Pomerons in the considered diagram. 




Figure 2: (a) Cylindrical diagram corresponding to the one-Pomeron exchange contribution to elastic 
pp scattering, and (b) the cut of this diagram which determines the contribution to the inelastic pp 
cross section. Quarks are shown by solid curves and string junction by dashed curves. 

For a nucleon target, the inclusive rapidity y (or Feynman-x xp) spectrum of a 
secondary hadron h has the form [1]: 

dfi X j7] dcT dfi ^ — ^ , h /■ \ I h / \ /\ 

— = -— = — =\^Wn-rn{x)+WD-rD{x), 8 

dy ainei dxp dy 

where the functions 0^(a;) determine the contribution of the diagram with n cut 
Pomerons and Wn is the relative weight of this diagram. The last term wd ■ 4>%{x) 
accounts for the contribution of diffraction dissociation processes. 

For pp collisions 



<l>pp{x) = f!;g{x+,n)- fl^{x-,n) + f!^{x+,n)- f!;g{x.,n) + 

+ 2-{n-l)-f^{x+,n)-f^{x^,n), (9) 



x± 



1 
2 



Arrij, /s + x"^ ± X 



(10) 



where fgg, fg, and fs correspond to the contributions of diquarks, valence quarks, and 
sea quarks, respectively. 
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These functions are determined by the convolution of the diquark and quark dis- 
tributions (that are normahzed to one) with the fragmentation functions, e.g. for the 
quark one can write: 

1 

fg{x+,n) = J Ug{xi,n) ■ G'^{x+/xi)dxi . (11) 

x+ 

At very high energies both x+ and X- are neghgibly small in the midrapidity region, 
and so all fragmentation functions, which are usually written [13] as C^i^z) = a/i-(l— 2;)^, 
become constants and equal for a particle and its antiparticle: 

G'J(a;+/xi) = ah . (12) 

This leads, in agreement with [23], to 

^ = 9H- {s/sor^'^-' ~ al . is/soT^^'^-' , (13) 
dy 

that corresponds to the only one-Pomeron exchange diagram (AGK theorem [12]) at 
asymptotically high energy. The values of the Pomeron parameters presented in Eq. (7) 
are used in the QGSM numerical calculations. 

The values dn/dy in the central region can be obtained under different conditions. 
Sometimes they are presented for all inelascic interactions, while in other cases they 
correspond to the events without single diffraction (NSD), or to events called INEL> 0, 
in which as minimum one charged particle should be detected in the kinematical window 
|?7| > 1 [25j. When considering the experimental value dn/dy = N particles /N events one 
has to keep in mind that though for both NSD and INEL> triggers the value of 
N particles in midrapidity region at LHC energies is constant, the number of events 
N events significantly changes in these three cases, leading to different values of dn/dy. 

As an example, we present in Table 2 three values of dn/dy measured by the ALICE 
Collaboration [21 [251 126] at ^/s = 900 GeV. 





All inelastic [24J 


NSD 


24 




INEL> 


'251 [26] 


dn/dy 


3.02 ±0.01°:°^ 


3.58±0.01°:i2 


3.81 ±0.01^:^^ 



Table 2. The inclusive densities of charged secondaries measured by ALICE Collaboration in all 
inelastic events, in events without single diffraction (NSD), and in events with as minimum one charged 
particle in the kinematical window \'ri\ > 1 (INEL> 0), in pp collisions at ^ = 900 GeV [24l [25l [26] . 
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Figure 3: The QGSM predictions for the pseudorapidity distributions of ah charged secondaries 
produced in inelastic pp and pp cohisions at different energies. 



The QGSM predictions for the dn/drj distributions of all charged secondaries pro- 
duced in inelastic pp and pp collisions at different energies are shown in Fig. 3. The 
experimental data are taken from ref . |23| . 

The QGSM allows one to calculate the inclusive spectra of different secondaries. 
The comparison of our predictions [8] with new data by the ALICE Collaboration [28] 
is presented in Table 3. 



Particle 


QGSM [8] 


ALICE Collaboration 


[28 




p 
A 

n+ 


1.68 

1.66 

0.17 

0.16 

0.10 

0.05 

0.005 

0.0004 


1.493 ±0.004 ±0.074 
1.485 ±0.004 ±0.074 
0.184 ±0.004 ±0.015 
0.183 ±0.004 ±0.015 
0.077 ±0.002 ±0.006 
0.08 
0.009 
0.0008 



Table 3. The QGSM predictions [8] for the midrapidity yields dn/dy {\y\ < 0.5) of different secon- 
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daries at energy y/s = 900 GcV together with the corresponding experimental data by the ALICE 
Cohaboration }28) . 

The inclusive densities of pions are slightly overestimated. The predicted |8J value 
of the ratio of H~/A midrapidity yields (~ 0.10) is in agreement with the experimental 
result [21] S~/A ~ 0.11 ± 0.005. The agreement of the order of 10—15% for all cases 
can be considered reasonable one. 

We predict [S] for all inelastic interactions one increase of the midrapidity yields 
for 71 and K mesons of about 1.4 times in the energy region ^/s = 900 GeV— 7 TeV 
(smaller than the increase in the data [30] for equal to 1.69 ± 0.01 ± 0.06 in NSD 
events). For antibaryons we predict [8] increases from ~ 1.6 for p and A to ~ 2.0 for 

n. 

4 Baryon/antibaryon asymmetry in QGSM 

In the string models baryons are considered as configurations consisting of three con- 
nected strings (related to three valence quarks) called string junction (SJ) [311 |32l [331 
[31], as it is shown in Fig. 4. Such a baryon structure is supported by lattice calculations 




Fi gure 4: The composite structure of a baryon in string models. Quarks are shown by open points 
and SJ by black point. 

This picture leads to some general phenomenological predictions. In particular, it 
opens room for exotic states, such as the multiquark bound states, 4-quark mesons, 
and pentaquarks [331 (Ml [37] . In the case of inclusive reactions, the baryon number 
transfer to large rapidity distances in hadron-nucleon reactions can be explained [5] by 
SJ diffusion. 
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The production of a baryon-antibaryon pair in the central region usually occurs via 
SJ-SJ pair production. Then, the SJ, that has upper color indices, combines with sea 
quarks, while the antiSJ [SJ), that has lower indices, combines with sea antiquarks, 
resulting, respectively, into the B and the B of a. BB pair [33] 138] , as it is shown in 
Figure 5a. 

d B 



c 




b c d 



Fi gure 5: QGSM diagrams describing secondary baryon B production by diquark d. (a) Central 
production of BB pair via production of new SJ-SJ pair. Single production of B in the processes of 
diquark fragmentation: (b) initial SJ together with two valence quarks and one sea quark, (c) initial 
SJ together with one valence quark and two sea quarks, and (d) initial SJ together with three sea 
quarks. Quarks are shown by solid curves and SJ by dashed curves. 

In processes with incident baryons, e.g. in pp collisions, another possibility to 
produce a secondary baryon in the central region exists. This second possibility is 
connected with the diffusion in rapidity space of a SJ existing in the initial state 
and it leads to significant differences in the yields of baryons and antibaryons in the 
midrapidity region even at rather high energies [5]. Probably, the most important 
experimental fact in favour for this process is the rather large asymmetry in Q and fl 
baryon production in high energy 7T~p interactions |39] . 

The quantitative theoretical description of the baryon number transfer via S J mech- 
anism was suggested in the 90 's and used to predict |10] the p/p asymmetry at HERA 
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energies. The quantitative description of the baryon number transfer due to SJ diffu- 
sion in rapidity space was first obtained in and then in papers P EJ [HI SH 1121 113] • 

To obtain the net baryon charge we consider three different possibihties according 
to ref. in]. The first one is the fragmentation of the incoming diquark giving rise to 
a leading baryon (Figure 5b). A second possibihty is to produce a leading meson in 
the first break-up of the string and one baryon in a subsequent break-up (Figure 5c). 
In these two first cases the baryon number transfer is possible only for short distances 
in rapidity. In the third case, shown in Figure 5d, both initial valence quarks in the 
diquark recombine with sea antiquarks into mesons M, while a secondary baryon is 
formed by the SJ together with three sea quarks. 

The fragmentation functions for the secondary baryon B production corresponding 
to the three processes shown in Figs. 5b-d can be written as follows (see [5J for more 
details): 

(14) 

z^-{l-z), (15) 
fs ■ z'-'^'-' ■ (1 - zf , (16) 







= ttN 








= ttN 








= aN 


■ e ■ v_ 



for Figs. 5b, 5c, and 5d, respectively, and where aj^ is the normalization parameter 
that was determined from the experimental data at fixed target energies, and f^, f^, 
w£ are the relative probabilities for different baryons production that can be found by 
simple quark combinatorics HS] . 

The contribution shown in Figure 5d is essential if the intercept of the S J exchange 
Regge-trajectory, asj, is large enough. The contribution of the graph in Figure 5d 
is weighted in QGSM by the coefficient e which determines the small probability for 
such a baryon number transfer to occur. Let's finally note that this process can be 
very naturally realized in the quark combinatorial approach [33] through the specific 
probabilities of a valence quark recombination (fusion) with sea quarks and antiquarks. 

At high energies the SJ contribution to the inclusive cross section of secondary 
baryon production at large rapidity distance Ay from the incident nucleon can be 
estimated as 

^ as ■ £ ■ e"(l-"^-)•^^ (17) 



cr dy 
where as = cln ■ vf^. 

The data by the ALICE Collaboration [36] for p/p ratios in pp collisions at LHC 
energies y/s = 900 GeV and 7 TeV are presented in Table 4, together with QGSM the 
corresponding calculations with different values of asj- These data seem to favour the 
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value asj = 0.5, what could mean that the Odderon contribution is not seen [7j in this 
process. 



SJ exchange 


^ = 900 GeV 


v/i = 7 TeV 


asj = 0.9 


0.89 


0.95 


asj = 0.5 


0.95 


0.99 


e = 


0.98 


1. 


ALICE [46J 


0.957 ±0.006 ±0.014 


0.991 ±0.005 ±0.014 



Table 4. The QGSM predictions for p/p in pp collisions at LHC energies and the data by the ALICE 
Collaboration ^46^. The value e = corresponds to the case with not C-negative exchange. 

The LHCb Collaboration measured the ratios of A to A in the rapidity interval 
2 < y < 4 at a/s = 900 GeV and 7 TeV |17]. We compare these results with the 
QGSM calculations in Fig. 6. 




Figure 6: The QGSM predictions (solid curves for asj = 0.9 and dashed curves for asj = 0.5) for 
the ratios of the spectra of secondary A to A as the functions of their rapitities at energies y/s = 900 
GeV (left) and i/s = 7 TeV (right), together with the experimental data by the LHCb Collaboration 

One can see that here the QGSM calculation with the value asj = 0.9 is in a slightly 
better agreement with the data. It is also necessary to note that the predictions of 
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PYTHIA 6 MC generator are in disagreement with the experimental data at i/i = 900 
GeV. 

We predict practically equal B /B ratios for baryons with different strangeness con- 
tent. 

5 Conclusion 

The first experimental data obtained at LHC are in general agreement with the calcu- 
lations provided by QGSM in the framework of Regge theory with the same values of 
parameters that were determined at lower energies (mainly for the description of the 
data of fixed target experiments). 

We neglect the possibility of interactions between Pomerons (so-called enhancement 
diagrams) in the calculations of integrated cross sections and inclusive densities. Such 
interactions are very important in the cases of heavy ion and nucleon-nucleus 
[19] interactions at RHIC energies, and their contribution should increase with energy. 
However, we estimate ^9] that the contributions of these enhanced diagrams to the 
inclusive density of secondaries produced in pp collisions at LHC energies is not large 
enough to be significant. 
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